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ABSTRACT 

The computer program CARES/LIFE calculates the time- 
dependent reliability of monolithic ceramic components subjected 
to therm omechanic al and/or proof test loading. This program is 
an extension of the CARES (Ceramics Analysis and Reliability 
Evaluation of Structures) computer program. CARES/LIFE 
accounts for the phenomenon of subcritical crack growth (SCG) 
by utilizing the power law, Paris law, or Walker equation. The 
two-parameter Weibull cumulative distribution function is used 
to characterize the variation in component strength. The effects 
of multiaxial stresses are modeled using either the principle of 
independent action (PIA), the Weibull normal stress averaging 
method (NSA), or the Batdorf theory. Inert strength and fatigue 
parameters are estimated from rupture strength data of naturally 
flawed specimens loaded in static, dynamic, or cyclic fatigue. 
Application of this design methodology is demonstrated using 
experimental data from alumina bar and disk flexure specimens 
which exhibit SCG when exposed to water. 


NOMENCLATURE 

A surface area; material-environmental fatigue constant 

a crack half length 

B subcritical crack growth constant 

C Shetty’s constant in mixed-mode fracture criterion 

g g-factor 

H step function 


i ranking of ordered fracture data in statistical analysis 

K stress intensity factor 

k crack density coefficient 

m Weibull modulus, or shape parameter 

N material-environmental fatigue constant 

n number of cycles 

P f cumulative failure probability 

Q cyclic fatigue parameter 

R ratio of minimum to maximum effective 

stress in a load cycle 
T period of one cycle 

t time 

tjj time-dependent scale parameter 

V volume; crack velocity 

x,y,z Cartesian coordinate directions 

Y crack geometry factor 

a angle between a n and the stress a 1 

P angle between a n projection and the stress c 2 in 

plane perpendicular to Oj 
A increment 

T| crack density function 

n 3.1416 

a applied stress distribution 

a 0 Weibull scale parameter 

a, a 2 ,o 3 tensor stress components; principal stresses 
(o, ;> o 2 > ct 3 ) 
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X 

shear stress acting on oblique plane whose normal is 
determined by angles a and Ji 

'F 

spatial location (x,y,z) and orientation (a,p) in a 
component 

n 

solid angle in three-dimensional principal stress 
space for which a e > C CT 

CO 

angle in two-dimensional principal stress space for 
which a e > a cr 

Subscripts: 

B 

Batdorf 

c 

cyclic; critical 

ch 

characteristic 

cr 

critical 

d 

dynamic fatigue 

e,ef 

effective 

eq 

equivalent 

f 

failure; fracture 

I 

crack opening mode 

n 

crack sliding mode 

m 

crack tearing mode 

i 

i’th 

max 

maximum 

min 

minimum 

n 

normal; normal stress averaging 

S 

surface 

T 

transformed 

u 

uniaxial 

V 

volume 

w 

Weibull 

0 

characteristic 


Superscripts: 

- modified parameter 

normalized quantity 


INTRODUCTION 

Advanced ceramic components designed for gasoline, diesel, 
and turbine heat engines are leading to lower engine emissions, 
higher fuel efficiency, and more compact designs due to their 
low density and ability to retain strength at high temperatures. 
Ceramic materials are also used for wear parts (nozzles, valves, 
seals, etc.), cutting tools, grinding wheels, bearings, coatings, 
electronics, and human prostheses. Among the many require- 
ments for the successful application of advanced ceramics are the 
proper characterization of material properties and the use of a 
mature and validated brittle material design methodology. 

Ceramics are brittle and the lack of ductility leads to low strain 
tolerance, low fracture toughness, and large variations in 
observed fracture strength. The material as processed has 
numerous inherent randomly distributed flaws. The observed 
scatter in fracture strength is caused by the variable severity of 
these flaws. The ability of a ceramic component to sustain a load 
also degrades over time. This is due to a variety of effects such 
as oxidation, creep, stress corrosion, and cyclic fatigue. Stress 
corrosion and cyclic fatigue result in a phenomenon called 
subcritical crack growth (SCG). SCG initiates at an existing flaw 
and continues until a critical length is reached, causing cata- 


strophic crack propagation. The SCG failure mechanism is a 
load-induced phenomenon over time. It can also be a function of 
chemical reaction, environment, debris wedging near the crack 
tip, and deterioration of -bridging ligaments. Because of this 
complexity, models that have been developed tend to be 
semi-empirical and approximate the behavior of subcritical crack 
growth phenomenologically. 

The objective of this paper is to present a description of the 
integrated design computer program, CARES/LIFE^ (Ceramics 
Analysis and Reliability Evaluation of Structures LIFE 
prediction program). The theory and concepts presented in this 
paper reflect the capabilities of the CARES/LIFE program for 
time-dependent probabilistic design. To determine the validity of 
the design methodology utilized in this software, time -dependent 
reliability predictions from CARES/LIFE are compared with 
experimental data^ from uniaxially and biaxially loaded 
alumina flexure bars and disks which are known to exhibit slow 
crack growth in water. 

PROGRAM CAPABILITY AND DESCRIPTION 

Probabilistic component design involves predicting the 
probability of failure for a thermomechanically loaded 
component from specimen rupture data. Typically these 
experiments are performed using many simple geometry flexural 
or tensile test specimens. A static, dynamic, or cyclic load is 
applied to each specimen until fracture. Statistical strength and 
SCG (fatigue) parameters are then determined from these data. 
Using these parameters and die results obtained from a finite 
element analysis, the time-dependent reliability for a complex 
component geometry and loading is then predicted. Appropriate 
design changes are made until an acceptable probability of 
failure has been reached. This design methodology combines the 
statistical nature of strength-controlling flaws with the mechanics 
of crack growth to allow for multiaxial stress states, concurrent 
(simultaneously occurring) flaw populations, and scaling effects. 
These issues are addressed within the CARES/LIFE program. 

CARES/LIFE predicts the probability of failure of a 
monolithic ceramic component as a function of service time. It 
assesses the risk that the component will fracture prematurely as 
a result of subcritical crack growth. The effect of proof testing 
components prior to service is also considered. CARES/LIFE is 
coupled to commercially available finite element programs such 
as ANSYS, ABAQUS, MSC/NASTRAN, and COSMOS/M. 
CARES/LIFE is an extension of the CARES (3,4) program. It 
retains all of the capabilities of the previous CARES code, 
which include fast-fracture component reliability evaluation and 
Weibull parameter estimation from inert strength (without SCG 
contributing to failure) specimen data. CARES/LIFE can 
estimate parameters that characterize SCG from specimen data 
as well. 

Finite element heat transfer and linear-elastic stress analyses 
are used to determine the component’s temperature and stress 
distributions. The reliability at each element is calculated 
assuming that randomly distributed volume flaws and/or surface 
flaws control the failure response. The probability of survival for 
each element is assumed to be a mutually exclusive event, and 
the overall component reliability is then the product of all the 
element survival probabilities. CARES/LIFE generates a data file 
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containing element risk-of- rupture intensities (a local measure of 
reliability) for graphical rendering of the structure’s critical 
regions. 

CARES/LIFE describes the probabilistic nature of material 
strength, using the Weibull cumulative distribution function. The 
effect of multiaxial stresses on reliability is predicted using the 
principle of independent action (PIA)/ the Weibull normal 
stress averaging method (NSA),^ or the Batdorf theory/ 8 * The 
Batdorf theory combines the weakest link theory and linear 
elastic fracture mechanics (LEFM). Conventional fracture 
mechanics analysis requires that both the size of the critical 
crack and its orientation relative to the applied loads determine 
the fracture stress. The Batdorf theory includes the calculation of 
the combined probability of the critical flaw being within a 
certain size range and being located and oriented so that it may 
cause fracture. A user-selected flaw geometry and a mixed-mode 
fracture criterion are required to model volume- or surface- 
strength-limiting defects. Mixed-mode fracture refers to the 
ability of a crack to grow under the combined actions of a 
normal load (opening mode) and shear Load (sliding and tearing 
modes) on the crack face. CARES/LIFE includes the total strain 
energy release rate fracture criterion, which assumes a crack will 
extend in its own plane (coplanar)/ 9 ^ Out-of-plane crack 
extension criteria are approximated by a simple semiempirical 
equation/ 10 ’ 11 ^ Available flaw geometries include the Griffith 
crack, penny-shaped crack, semicircular crack, and notched 
crack. The Batdorf theory is equivalent to the probabilistic 
multiaxial theories proposed by Evans^ 2 ^ and Matsuo/^ 

Subcritical crack growth is difficult to model, because it is a 
complex phenomenon often involving a combination of failure 
mechanisms. Existing models usually involve empirically derived 
crack propagation laws that describe the crack growth in terms 
of the stress intensity factor at the crack tip plus additional 
parameters obtained from experimental data. 

In CARES/LIFE, the relations describing subcritical crack 
growth are directly incorporated into the PIA, NSA, and Batdorf 
theories. Subcritical crack growth is modeled with the power 
l aw (l4,15) Paris law/ 16 -* and the Walker law* 17,18 ^ for static 
and constant-amplitude cyclic loading. These laws use 
experimentally determined parameters which are material- and 
environment-sensitive. The power law is used to model stress 
corrosion cracking in materials such as glasses and alumina 
exposed to H 2 0. Elevated-temperature slow crack growth of 
silicon nitrides, silicon carbides, and alumina also follows power 
law behavior. 

Some polycrystalline ceramics are prone to strength 
degradation due to mechanical damage induced by cyclic 
loading. The Paris and Walker laws have been suggested as 
models to account for this behavior/ 18 ^ The Walker equation is 
functionally similar to the Paris equation with additional terms 
to account for the effect of the R-ratio (minimum cycle stress to 
maximum cycle stress) on lifetime. 

CARES/LIFE is capable of predicting the change in a 
surviving component’s reliability after proof testing is performed. 
Proof testing is the loading of all components prior to service to 
eliminate those which may fail prematurely. The components 


that survive the proof test will have a lower (attenuated) risk of 
failure in service. In CARES/LIFE the attenuated failure 
probability is calculated using the PIA, the Weibull normal stress 
averaging, and the Batdorf theories. The Batdorf model is used 
to calculate the attenuated failure probability when the proof test 
load and the service load are not in line or have different 
multiaxial stress states. This feature is useful when the proof test 
does not identically simulate the actual service conditions on the 
component The durations of the proof test and the service load 
are also considered in the analysis. 

Predicted lifetime reliability of structural ceramic components 
depends on Weibull and fatigue parameters estimated from 
widely used tests involving flexural or tensile specimens. 
CARES/LIFE estimates fatigue parameters from naturally flawed 
specimens ruptured under static, cyclic, or dynamic (constant 
stress rate) loading. Fatigue and Weibull parameters are 
calculated from rupture data of three-point and four-point flexure 
bars, as well as tensile specimens. For other specimen 
geometries, a finite element model of the specimen is also 
required when estimating these parameters. 


THEORY 

Time-dependent reliability is based on the mode I equivalent 
stress distribution transformed to its equivalent stress distribution 
at time t=0. Investigations of mode I crack extension (19) have 
resulted in the following relationship for the equivalent mode I 
stress intensity factor 


KfeqOF.t) = OfcqCF.t) Y / aflM) (1) 

where a^OFtt) is the equivalent mode I stress on the crack, Y 
is a function of crack geometry, aCP,t) is the appropriate crack 
length, and ¥ represents a location (x,y,z) within the body and 
the orientation (a,P) of the crack. In some models such as the 
Weibull and PIA, ¥ represents a location only. Y is a function 
of crack geometry; however, herein it is assumed constant with 
subcritical crack growth. Crack growth as a function of the 
equivalent mode I stress intensity factor is assumed to follow a 
power law relationship 

da(y,t) = A Kj^CP.t) (2) 

dt ^ 


where A and N are material/environmental constants. The 
transformation of the equivalent stress distribution at the time of 
failure, t=tf, to its critical effective stress distribution at time t=0 
is expressed as^ 20 ’ 21 ^ 




jf'o^OP.O* 

B ^ 




CP.t f ) 


(3) 
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where 


AY 2 Kjc -2 (N-2) 

is a raaterial/environmental fatigue parameter, K lc is the critical 
stress intensity factor, and a^OF^) is the equivalent stress 
distribution in the component at time t=tf. The dimensionless 
fatigue parameter N is independent of fracture criterion. B is 
adjusted to satisfy die requirement that for a uniaxial stress state, 
all models produce the same probability of failure. The 
parameter B has units of stress 2 x time. 

Volume Raw Analysis 

The probability of failure for a ceramic component using the 
Batdorf model (8,9,22) for volume flaws is 


PfV = 1 “ ex P r 


-f 

°T a 

dV 

J 

V 

{ 4k dc w CT ! 



(4) 


where V is the volume, T| v is the crack density function, a c 
is the maximum value of 0 f° r values of Y, and is tKe 
area of a solid angle projected onto a unit radius sphere in 
principal stress space containing all crack orientations for which 
the effective stress is greater than or equal to the critical mode 
I strength, a ^ The crack density distribution is a function of the 
critical effective stress distribution. For volume flaw analysis, the 
crack density function is expressed as 


Hv^crC**)) = k B V 


(5) 


Pfv(tf) = 1 -e *P 


■■if JiT-T °S.o cP ' t f ) sina da dp dV 

(7) 


Fracture criteria and crack shapes available for time-dependent 
analysis are identical to those used for fast fracture analysis in 
CARES.®'** These fracture criteria include Weibull normal 
stress averaging (a shear-insensitive case of the Batdorf theory), 
the total coplanar strain energy release rate, and the noncoplanar 
crack extension (Shetty) criterion. 

For a stressed component, the probability of failure for volume 
flaw analysis is calculated from equation (7). The finite element 
method enables discretization of the component into incremental 
volume elements. CARES/LIFE evaluates the reliability at the 
Gaussian integration points of the element or, optionally, at the 
element centroid. Subelement volume is defined as the 
contribution of the integration point to the element volume in the 
course of the numerical integration procedure. The volume of 
each subelement (corresponding to a Gauss integration point) is 
calculated using shape functions inherent to the element type^. 
Assuming that the probability of survival for each element is a 
mutually exclusive event, the overall component reliability is 
then the product of all the calculated element (or subelement) 
survival probabilities. 


Surface Raw Analysis 

The probability of failure for a ceramic component using the 
Batdorf model^’ 9,22) for surface flaws is 


p fs = l- ex P 



CO <fr|s( a cr) 

da cr 


do„ 


dAf 


( 8 ) 


where k BV and m v are material constants. The solid angle is 
expressed as 

Q = So* So H(o feq.0 ,o cr) sina da d P (6) 

where 


where A is die surface area, r\ s is the crack density function, 
a e is the maximum value of Gj^ q for all values of *F, and to 
is file arc length of an angle oc projected onto a unit radius semi- 
circle in principal stress space containing all of the crack 
orientations for which the effective stress is greater than or equal 
to tire critical stress. Analogous to the argument for volume 
flaws, equation (8) can be reformulated, yielding 1 ' 22 ' 1 


H^Ieq^cr) = 


1 Oleq,0^cr 
0 <*leq,0< a cr 


p fs( l f) = 1 " CX P 



IDe 

J leq,0 


0F,t f ) da dA 


(9) 


and a and P are the radial and azimuthal angles, respectively, on 
the unit radius sphere. The transformed equivalent stress 
is dependent on the appropriate fracture criterion, crack shape, 
and time to failure, if. Equation (4) can be simplified by 
performing the integration of yielding the timendependent 
probability of failure for volume flaw analysis 


The transformed equivalent stress G Ie q 0 is dependent on the 
appropriate fracture criterion, crack shape, and time to failure, 
tf. The fracture criteria and crack shapes available for time- 
dependent analysis are identical to those used for fast fracture 
analysis. These fracture criteria include Weibull normal stress 
averaging (a shear-insensitive case of the Batdorf theory), die 
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total coplanar strain energy release rate, and the noncoplanar 
crack extension (Shetty) criterion. 

The Finite element method enables discretization of the surface 
of the component into incremental area elements. CARES/LIFE 
evaluates the failure probability at the Gaussian integration 
points of shell elements or, optionally, at the element centroids. 
The area of each subelement (corresponding to a Gaussian 
integration point) is calculated using shape functions inherent to 
the element type^. Assuming that the probability of survival for 
each element is a mutually exclusive event, the overall 
component reliability is then the product of all the calculated 
element (or subelement) survival probabilities. 


Static Fatigue 

Static fatigue is defined as the application of a nonvarying load 
over time. For this case the mode I equivalent stress, OfcqCP.0. 
is independent of time and is thus denoted by o^OF)- 
Integrating equation (3) with respect to time yields 


<W’ P - t f) * OfcqW 


tf4qW 

B 


+ 1 


"Ur 


( 10 ) 


Dynamic Fatigue 

Dynamic fatigue is defined as the application of a constant 
stress rate <*0P) over a period of time, t Assuming the applied 
stress is zero at time t=0, then 


OfcqCF.t) = <W t 


(ID 


stress, q Tfr CF,t), at some location 'P. The equivalent static stress 
is related to the cyclic stress by 


<&,('*'> ti = <£,0^) dt = ‘i 



(13) 


= gW Ofcqc, 


m t! 


The CARES/LIFE program uses the maximum cyclic stress, 
a Ieqc ^OF), of the periodic load as a characteristic value to 

normalize the g-factor. For a periodic load over a time t lt the 
mode I static equivalent stress distribution is 


l 


(14) 

The use of g-factors for determining component life is an 
unconservative practice for materials prone to cyclic damage. 
The Walker equation/ 17 ^ which has traditionally been used in 
metals design, has been suggested as a model of fatigue damage 
for some ceramic materials/ 18) The Walker equation describes 
the crack growth increment per cycle, n, as 

da0y>n) = A K k Q (*P»n) AK^CP.n) (15) 
dn 


= °Ieqc Iiax W 


gwt f m 




B 


+ 1 


Substituting equation (11) into equation (3) results in an 
expression for effective stress at the time of failure 


where 

K kmM CP,n) = o IeqCmu 0P.n) Y Rm 




CP) = 


'1 <• <£>.■,) 

(N+l) B v f 


1 


(12) 


Cyclic Fatigue 

Cyclic fatigue is defined as the repeated application of a 
loading sequence. Analysis of the time-dependent probability of 
failure for a component subjected to various cyclic boundary 
load conditions is simplified by transforming that type of loading 
to an equivalent static load. The conversion satisfies the 
requirement that both systems will cause the same crack 
growth/ 23 ^ Implicit in this conversion is the validity of equation 
(2) for describing the crack growth. The probability of failure is 
obtained with respect to the equivalent static state. 

Evans (24) and Mencik (23) defmed g-factors, gOF), for various 
types of cyclic loading, that are used to convert the cyclic load 
pattern to an equivalent static load. For periodic loading, T is the 
time interval of one cycle, and cr^OF) is the equivalent static 
stress acting over the same time interval, t 1? as the applied cyclic 


and AK k CF,n) represents the range of the stress intensity factor 
over the load cycle. The subscripts max and min indicate the 
maximum and minimum cycle stress, respectively. The cyclic 
fatigue parameters A, N, and Q are experimentally determined. 
The Walker equation reduces to the Paris law (l6 ^ when N and Q 
are equal in value. The integration of Eq. (15) parallels that of 
Eq. (2), yielding the cyclic fatigue equivalent stress distribution 


Oleqc.O^f) 


where R0F,n) is the ratio of the minimum to maximum cyclic 
stress, n^ is the number of cycles to failure, and B is now 
expressed in units of stress 2 x cycle. The parameters B and N 
are determined from cyclic data. 
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1/(N S +1) 


Evaluation Of Fatigue Parameters from Inherently 
Flawed Specimens 

Lifetime reliability of structural ceramic components depends 
on the history of the loading, the component geometry, the 
distribution of pre-existing flaws, and the parameters that 
characterize subcritical crack growth. These crack growth 
parameters must be measured tinder conditions representative of 
the service environment When determining fatigue parameters 
from rupture data of naturally flawed specimens, the statistical 
effects of the flaw distribution must be considered along with the 
strength degradation effects of subcritical crack growth. In the 
following discussion, fatigue parameter estimation methods are 
described for surface flaw analysis using the power law 
formulation for constant stress rate loading (dynamic fatigue). 
Analogous formulations for volume flaws, static fatigue, and 
cyclic fatigue have also been developed/ 1 ^ 

For die uniaxial Weibull distribution the probability of failure 
is expressed as 


PfsCtf) = 1 - exp 


-kws J dA 

A 


(17) 


where 4* represents a location (x,y) and o 10 denotes the 
transformed uniaxial stress analogous to Ow 0 as defined in 
equation (13). The Weibull crack density coefficient is given by 


A d " 


Nc-2 

(N s + 1) B wS c oS 


V 


1-0.50 JJ 


l/m s 


where 


IHc 


ro c — 


N s -2 


( 21 ) 


The constants and Agf are obtained by equating risks of 
rupture. A cf is a modified effective area required for the time- 
dependent formulation. For the uniaxial Weibull distribution, the 
expression for the modified effective area is 


f 

'ojOF-tf)' 

J 

A 

G f 

l 1 ) 


dA 


( 22 ) 


where G { Q¥yif) denotes the maximum principal stress 
distribution. For multiaxially stressed components, the Batdorf 
technique is used to evaluate fatigue parameters. The analogous 
formulation for A cf is then 


k ws = 4 r < 18 > 

°oS 

The Weibull scale parameter, a oS , corresponds to the stress level 
where 63.21 percent of specimens with unit area would fail and 
has units of stress x area Vm *. CARES/LIFE normalizes the 
various fracture criteria to yield an identical probability of failure 
for the uniaxial stress state. This is achieved by adjusting the 
fatigue constant B. For the uniaxial Weibull model this adjusted 
value is denoted by B w and for the Batdorf model it is denoted 
by B b . From the dynamic fatigue equation (12), substituting B wS 
for B, N s for N, the uniaxial stress g x for G^, and rearranging 
equation (17) while assuming that 

o^(T,t f ) t f ^ i (19) 

(N s+ 1) B wS 


r 

rf = — J 


f 


a Ieq (T,t f ) 


Wn, 


l s^s 


da 


dA 


(23) 


where die normalized Batdorf crack density coefficient is 


k BS = k BS /k wS . Equation (21) is applicable except that B bs 
replaces B w c. The relationship between B B c and B wS for a 
uniaxial load* is established by equating the risk of rupture of die 
Batdorf model with that of the uniaxial Weibull model (1) 


wS 

BS 


r mcNc 

it jo, s s 0F,t f ) dA 


moNr 


2 k BS /X 7 S CP. £ f) da dA 


l^Tlc 


(24) 


the median behavior of the experimental dynamic fatigue data 
can be described by 


*0.5 


= 


l/(N s +l) 


( 20 ) 


where G^ is the median rupture stress of the specimen and d 
represents the stress rate at the location of maximum stress. The 
constant A d is 


As N s becomes large, equation (24) approaches unity. 

The terms A^ and N s in equation (21) are determined from 
experimental data. Taking the logarithm of equation (20) yields 

In <J f = In A d + L_ (net ( 25 > 

d N s + 1 

Linear regression analysis of the experimental data is used to 
solve equation (25). The median value method is based on least 
squares linear regression of median data points for various stress 


6 


Noel N. Nemeth 



rates. Another technique uses least squares linear regression on 
all the data points. A third option for estimating fatigue 
parameters is a modification to a method used by Jakus/ 2 ^ In 
this procedure, fatigue parameters are determined by minimizing 
the median deviation of the logarithm of the failure stress. The 
median deviation is the mean of the residuals, where the residual 
is defined as the absolute value of the difference between the 
logarithm of the failure stress and the logarithm of the median 
value. In CARES/LIFE this minimization is accomplished by 
maximizing ro<j (N<j + 1) estimated from the data versus the 
fatigue exponent 

To obtain A^ based on the median line of the distribution the 
following steps are taken. Experimental data at a sufficient 
number of discrete levels of applied stressing rate are 
transformed to equivalent failure times t Ti at a fixed stress rate 
fcj. (equating failure probabilities using equation (17)) 


l Ti 




r«. 



"VCNs+n 


(26) 


where the subscript T indicates a transformed value, the 
subscript i denotes each observed data number, and 
l fi = 1x1 CARES/LIFE the value of 6 X is taken as the 

lowest level of stressing rate in the data set With the data 
defined by a single Weibull distribution, parameter estimation is 
performed on the transformed data using 


Pfi 



r 


Ti 


V*Uj(N s +l) 



(27) 


where die characteristic time is 


l/CNj+l) 

(28) 


CARES/LIFE performs least squares or maximum likelihood 
Weibull parameter estimation as described in Pai*^ to solve 
equation (27) for m^ (Ng + 1) and tjQg. Substituting for 
fc in equation (20) and solving for the rupture stress in equation 
(27) corresponding to a 50% probability of failure, yields 


]1/(N S +1) 



1 T8S 


No-2 


(N s + 1) B oj 


T 


A C f “o 


This value of A d is used with the fatigue exponent N s estimated 
either with the least squares (using all experimental rupture 
stresses) or median deviation method. The fatigue constant B is 
obtained from equations (21) or (28). 


EXAMPLE 

This example demonstrates the use of CARES/LIFE to predict 
the time-dependent reliability of components under multiaxial 
loads. The data for this example is from experimental work 
performed by Chao and Shetty (2) on alumina disks and bars 
exposed to water at room temperature. Chao and Shetty 
investigated the relationship between stress state and time- 
dependent strength degradation, specifically to determine if 
strength degradation due to slow crack growth in biaxial flexure 
can be predicted from inert fracture stresses and dynamic fatigue 
assessed in simple uniaxial tests. Time-dependent reliability 
predictions for the alumina from CARES/LIFE are compared to 
the results obtained by Chao and Shetty. 

Details regarding specimen preparation and testing are given 
in reference 2. Two batches of an alumina ceramic were pur- 
chased from the same vendor. The first batch was in the form of 
plates (127x127x5 mm) and rods (50.8 mm diameter, 76.2 mm 
length). The plates and rods were made from the same powder 
lot with identical isostatic pressing and sintering conditions. The 
second batch of alumina, in the form of rods, was purchased 
subsequently to examine dynamic fatigue under biaxial stresses. 
This material had identical chemistry and preparation as the first 
batch. The measured properties of the first batch include a 
fracture toughness, K^, of 4.13 MPa • { m , Young’s modulus, 
E, of 297.2 GPa, and a Poisson’s ratio, v, of 0.23. 

Bar specimens were cut from the plate stock and disks were 
cut from the rod stock. Specimens were either tested in a dry 
nitrogen environment at 100 MPa/s to obtain inert (fast-fracture) 
strengths, or in deionized water at various stress rates to obtain 
the time-dependent fracture strengths. The specimens were 
carefully prepared to minimize machining damage or failure 
from edge flaws. The uniaxial specimens were nominally 
3x4x45 mm. These specimens were loaded in either four-point 
flexure with an outer support span of 40 mm and an inner load 
span of 20 mm, or in three-point flexure using a 40 mm support 
span. The disk specimens were nominally 3.175 mm in thickness 
and 50.8 mm in diameter. They were loaded under uniform 
pressure on one face and supported on the other face by 40 
freely rotating ball bearings spaced uniformly along a 49.53 mm 
diameter circle. Fractography showed that all specimens broke 
due to a single population of randomly oriented surface flaws. 

Table 1 lists the fast-fracture Weibull parameters, estimated 
using the maximum likelihood method, for the various specimen 
configurations and loads. The values shown correlate very well 
to those of reference 2 for the Weibull shape parameter and 
characteristic strength. The results for the 90% confidence 
intervals differ somewhat from reference 2 due to the methods 
of estimation. CARES/LIFE uses the technique from Thoman, 
Bain and Antle (27) , while reference 2 uses a bootstrap technique. 
The 90% confidence interval on the Weibull modulus 
significantly overlaps for both bar and disk tests. This is an 
indication that the strength response is controlled by the same 
flaw population. Confidence bands on the characteristic strength 
may only be compared for identical specimen loading and 
geometry due to the size effect. The confidence intervals 
obtained for the two disk tests using CARES/LIFE are shown to 
overlap. The large confidence interval for the batch 2 specimens 
is due to the small sample size of 7 specimens. In reference 2 
the discrepancies in average strength between batch 1 and batch 
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2 disks were attributed to material processing differences. This 
will be further discussed. 

Fatigue parameters were estimated from 4-point bend bar 
specimens loaded in dynamic fatigue in water. At least fifteen 
specimens each were tested at stress rates of 0.02, 0.1, 1.0, 10.0, 
and 100.0 MPa/s. Figure 1 shows a plot of the individual failure 
stresses of the 95 specimens tested. Superimposed on this figure 
are median lines and 5% and 95% confidence bounds on this 
data as estimated with the median deviation technique. Table 2 
lists the estimated fatigue parameters using the median value, 
least squares, and median deviation techniques. The median 
deviation technique shows best agreement with Chao and Shetty 
for this case. In reference 2, the power law is used in the 
following form 


da(t) _ v 
dt c 



\N 

/ 


(30) 


where V c is the critical crack velocity. The CARES/LIFE 
program is formulated using the fatigue constant B and, 
therefore, V c is not explicitly calculated. However, for 
comparison, V c can be computed from CARES/LIFE results 
using die fracture strength relation* 2 * 


G f = d 


2 k£ s (N s + 1) 

53 J* s 

V c Y Ns (N s - 2)<* 


' l-Nj/2 l-Ns/2! 
a i * a f ] 


1 

(Njny 


(31) 

where aj is initial crack size and a^ is the crack size at failure. 
The subscript S denotes surface-flaw-dependent properties. 
Assuming af » and rearranging the above expression yields 




G f = 


2 V (N s + 1) *i 
v c Y Ns (N s - 2) 


l-N s /2 


I l/CNj+l) 


> N+l 


(32) 


Equating equations (32) and (20) results in 


Nc 


l-Ns/2 


V c = 


2 K fc a (N s + 1) a ; 
A d Ns+1 Y Ns (N s - 2) 


(33) 


The initial crack length and fatigue constant A d are evaluated 
here for a failure probability of 50 percent The crack geometry 
factor Y for a semicircular crack is 1.366* 3 *. Crack length was 
determined from equation (1) using the strength determined from 
inert testing (Table 1) for a 50 percent probability of failure. 
Chao and Shetty* 2 * estimated that N s was 40.7 and V c was 9.1 
m/s for Y equal to 1.24. The differences between the various 
parameter estimates in CARES/LEFE and reference 2 are not 
considered significant The fatigue velocity V c is particularly 
sensitive to N s as shown in equation (33). Table 2 also lists the 


fatigue constants B w g and B B g. B bs is determined for a 
semicircular crack and_noncoplanar crack extension with a shear 
sensitivity constant of C = 0.82* 3 *. The differences between B wS 
and B bs are small since N s is relatively large. Finally, the 
Weibull modulus can be directly estimated from the fatigue data 
using the relation mg=mg/(Ng-2) when N$ and m§ are known. 
For example, using the median deviation method m s = 25.0, 
which is consistent with the results shown in Table 1. 

To obtain disk specimens for dynamic fatigue tests, a second 
batch of material (of identical dimensions as the first material 
batch) was secured by Chao and Shetty* 2 *. Seven of these 
specimens were broken under inert conditions and a total of 35 
specimens were broken at stressing rates of 0.02, 0.1, 10.0, and 
100.0 MPa/s. Table 1 lists the fast-fracture Weibull parameters 
of the (batch 2) disks and Table 3 gives the estimates for the 
fatigue parameters. Table 3 shows only B bs , since B wS is 
formulated only for the uniaxial stress state. Figure 2 shows a 
plot of the individual failure stresses as well as the median line 
and 5% and 95% confidence bounds estimated with the median 
deviation technique. One data point was flagged as an outlier (at 
100 MPa/s); however, it was not rejected and had little effect on 
the overall results. Batch 2 material showed a stronger than 
expected average strength relative to the batch 1 material Chao 
and Shetty* 2 * attributed this to a small increase of K k in the 
material. Although batch 2 was unexpectedly strong, the rate of 
strength degradation was similar to that of the batch 1 four-point 
flexure data. Reference 2 reported that the 90% confidence 
intervals of N s for both batches showed significant overlap. This 
was further confirmed in Table 4, which shows the ratio of the 
natural logarithm of the characteristic strength of the batch 2 
Hides to the natural logarithm of the characteristic strength of the 
batch 1 four-point flexure specimens at each stressing rate 
(maximum likelihood estimates). If the rate of change of strength 
degradation were stress-state dependent, then the strength ratio 
would systematically change with stressing rate. However, Table 
4 shows no such change. Therefore, the differences between 
fatigue exponents N s in Table 3 and Table 2 appear to be 
reflecting expected statistical variation. 

The effect of multiaxial stress states on the material is 
assessed by comparing the difference in inert strength between 
the uniaxially loaded four-point bend specimen and the biaxially 
loaded disk. Assuming that small crack-like imperfections 
control the failure, the material strength in multiaxial stress 
states can be correlated to the effects of mixed-mode loading on 
the individual cracks.**’^* Shetty* 1 ** developed a simple equation 
describing the ability of a crack to extend under the combined 
actions of a normal and shear loadjpn the crack face using an 
empirically determined parameter, C. For a semicircular crack 
this equation is* 3 * 


CT Ieq 


- 

' 

<*n + 

<*n + 3-301 fc? 

1 

[c) J 


(34) 


where a n and x are the normal and shear stresses, respectively, 
acting on the flaw plane. 

The failure response of the biaxial flexure specimens and the 
three-point flexure bars is predicted using the Weibull 
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parameters (Table 1) and fatigue parameters (Table 2; median 
deviation results) estimated from the four-point flexure bar 
rupture tests. The prediction for the three-point bars is compared 
to experi-mentai results to confirm the expected Weibull size 
effect The disk prediction is compared to experimental results 
in order to make assertions regarding the effect of multiaxial 
stress states on the fast-fracture and fatigue fracture of the 
material. 

For this analysis, the disk predictions are based on a finite 
element model of the disk as shown in Figure 3. Brick and 
wedge solid elements are used to model a 7.5° sector of the disk 
and appropriate boundary conditions are applied corresponding 
to a 49.53 mm diameter ring support and a pressure load on one 
face. A gradient is imposed on the nodal spacing such that die 
smallest elements have the highest tensile stresses. This is 
required to obtain a satisfactory convergence of the reliability 
solution and is independent of the mesh size needed to obtain 
accurate nodal stresses. To perform surface flaw reliability 
analysis quadrilateral and triangular shell elements are attached 
to the surface nodes on the disk tensile face in order to obtain 
surface stresses and areas. The shell elements are very thin and 
have membrane properties only. These elements are such that 
they contribute negligible stiffness to the model. Verification of 
the accuracy of the finite element model is obtained by com- 
parison to available closed form solutions. Chao and Shetty^ 
calibrated the applied pressure on the disk to the measured strain 
at the center of the disk. The strain calibration was used to 
correlate the fracture stress with the fracture pressure. The 
CARES/LIFE analysis uses these reported fracture stresses rather 
than the actual applied pressure on the disk. Reliability pre- 
dictions for the three-point flexure bar are obtained from closed 
form solutions of the effective area/ 1 ^ 

Probability of failure predictions are made for the disk and 
three-point bar in the fast-fracture condition and also for 
dynamic loading at a stressing rate of 1 MPa/s in water. The 
fracture strength distribution of dynamically (constant stress rate) 
loaded specimens is characterized by a Weibull distribution with 
a Weibull modulus, m^, of 


m dS 


N s + f 

Nc - 2 
V 5 ) 


m s 


(35) 


For this analysis m^ has a value of 25.58. Table 5 lists the 
estimated Weibull parameters obtained from bars and disks 
loaded in this condition. In all cases the estimated Weibull 
modulus is somewhat higher than predicted; however, the pre- 
dicted value was within the 90% confidence bounds. Table 5 
shows that the Weibull modulus did increase from the inert 
condition, as expected. If the Weibull modulus value of m s =25.0 
obtained directly from the four-point fatigue data were con- 
sidered, then would be 26.9, which is in better agreement 
with Table 5. 

The CARES/LIFE program is designed to predict reliability for 
static fatigue or constant-amplitude cyclic loading. To predict 
reliability for dynamic loading, an equivalent static loading time 
is computed using the approach given in equation (16). In this 
case, the g-f actor is equal to 1/(N S +1) multiplied by the time to 
failure. Figures 4 and 5 are the resulting Weibull plots of these 


predictions for various levels of shear sensitivity parameter C. 
In Figure 4 a value of C of 0.82 is used which corresponds to 
an approximation of the maximum tangential stress mixed-mode 
fracture criterion. The dashed line in the figure denotes the 
results from reference 2. The^ small differences indicated are 
mainly due to the different C and crack geometry factor, Y, 
used and are not due to the different values used for the fatigue 
parameters. This figure also shows a good correlation between 
predicted and experimental results for the three-point flexure 
bars, which confirms the size effect expected in fast-fracture and 
fatigue. The CARES/LIFE prediction lines for the three-point 
flexure bars superimpose on the three-point bar prediction lines 
of reference 2. _ 

Figure 5 shows predictions for C values of 0.82, 0.90, and 
»1 (solid, long dashed, and short dashed line, respectively). 
Since in this example the CARES/LIFE fracture predictions are 
normalized to the (uniaxial stress state) four-point flexure data, 
then the choice of a fracture criterion onlyjsffects the predictions 
for the (biaxial stress state) disks. The C values of 0.82 and 
much greater than one represent the extreme bounds on the 
expected mixed-mode shear sensitivity of the flaws in the 
material. The value of C»1.0 represents a shear-insensitive 
fracture criterion, while C=0.82 is_highly shear sensitive. The 
0.90 value represents a choice of C_that best fits to both inert 
and fatigue disk data. Optimizing C in this manner was not 
considered in reference 2. From this figure stress state effects are 
adequately accounted for, both in fast-fracture and subcntical 
crack growth conditions. 


CONCLUSION 

The use of structural ceramics for high-temperature 
applications depends on the strength, toughness, and reliability 
of these materials. Ceramic components can be designed for 
service if the factors that cause material failure are accounted 
for. This design methodology must combine the statistical nature 
of strength controlling flaws with fracture mechanics to allow 
for multiaxial stress states, concurrent flaw populations, and 
subcritical crack growth. This has been accomplished with the 
CARES/LIFE public domain computer program for predicting 
the time-dependent reliability of monolithic structural ceramic 
components. An example has been given to illustrate the use of 
this design methodology in CARES/LIFE for predicting the 
effects of component size, stress distribution, stress state, and 
subcritical crack growth on the lifetime reliability. Potential 
enhancements to the code include the capability for transient 
analysis, three-parameter Weibull statistics, creep and oxidation 
modeling, flaw anisotropy, threshold stress behavior, and 
parameter regression for multiple specimen sizes. 
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TABLE 1-WeibuII parameters estimated from inert data. 


Specimen 

Weibull 

modulus, 

m s 

90% confidence 
bounds on 

m s 

Characteristic 

strength, 

a es 

MPa 

90% confidence 
bounds on 

°es 

MPa 

Scale 

parameter, 

°oS 

MPa • mm 2/m * 

3- point 

4- point 
Disk 

Disk (batch 2) 

25.43 

23.76 

22.25 

28.98 

20.47, 29.91 
19.13, 27.95 
17.12, 26.81 
13.28, 40.88 

385.9 

353.4 
338.7 

351.4 

382.0, 389.9 
349.6, 357.3 

334.1, 343.5 
341.5, 362.2 

414.7 

425.8 
436.4 
423.2 


Estimation 

method 


Median Value 
Least Squares 
Med Deviation 
Reference 2 


TABLE 2-Fatigue parameters for 4-point bend bar^ 


Fatigue 

constant. 

Fatigue 

exponent, 

N s 

Fatigue 

constant, 

®wS 

MPa 2 • s 

Fatigue 

constant, 

®BS 

MPa 2 • s 

Crack 

velocity, 

V c 

m • s' 1 

2.339 x 10 2 

36.04 

4.590 x 10-' 

4.579 x 10-' 

1.162 

2.336 x 10 2 

40.84 

5.631 x 10* 

5.617 x 10* 

8.176 

2.337 x 10 2 

41.23 

4.783 x 10* 

4.771 x 10* 

9.373 

Mll 

40.7 

— 

— 

9.1 


TABLE 3— Fatigue parameters for disk (batch 2) specimens. 


Estimation 

method 

Fatigue 

constant, 

A d 

Fatigue 

exponent. 

Ns 

Fatigue 

constant, 

®BS 

MPa 2 * s 

Crack 

velocity, 

V e 

m • s" 1 

Median Value 

2.293 x 10 2 

35.68 

2.675 x 10* 

1.894 

Least Squares 

2.313 x 10 2 

41.79 

3.959 x 10* 

16.22 

Med Deviation 

2.304 x 10 2 

36.23 

2.649 x 10' 1 

1.988 

A A 

Reference 2 

— 

36.6 

— 

2.4 
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TABLE 4— Comparison of 4-point and disk (batch 2) 
characteristic strengths 


Stressing 

rate, 

MPa • s' 1 

4-point 

strength, 

s es 

MPa 

Disk 

strength, 

°0S 

MPa 

Strength 

ratio, 

to a 6S„ 


0.020 

2.145 x 10 2 

2.116 x 10 2 

0.9975 

0.100 

2.250 x 10 2 

2.195 x 10 2 

0.9954 

10.00 

2.531 x 10 2 

2.492 x 10 2 

0.9972 

100.0 

2.628 x 10 2 

2.620 x 10 2 

0.9995 

Inert 

3.534 x 10 2 

3.514 x 10 2 

0.9990 


TABLE 5-Weibull parameters estimated in water (1 MPa/sV 


Specimen 

Weibull 

modulus, 

m dS 

90% confidence 
bounds on 

m dS 

Characteristic 

strength, 

°es 

MPa 

90% confidence 
bounds on 

0 es 

MPa 

3-point 

27.83 

21.20, 33.70 

255.6 

252.7, 258.5 

4-point 

26.70 

20.54, 32.17 

236.3 

233.6, 239.0 

Disk 

32.66 

24.87, 39.53 

215.9 

213.8, 218.1 
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Figure 1 .—Dynamic fatigue of alumina four-point flexure bars in water. 
Median regression line (solid), and 5% and 95% regression lines 
(dashed) estimated with the median deviation technique. 



Figure 2. — Dynamic fatigue of alumina pressure-on-disk flexure speci- 
mens in water. Median regression line (solid), 5% and 95% regression 
fines (dashed) estimated with the median deviation technique. 



Figure 3. — Finite element model of a 7.5° sector of the disk. 
Arrows indicate pressure load on disk. 



Figure 4. — Weibull plot showing fast-fracture three-point and disk 
specimen strengths measured in an inert environment, as well as 
three-point and disk specimen fracture strengths dynamically 
loaded at a rate of 1 MPa/s. Solid lines are corresponding 
CARES/LIFE predictions for C = 0.82. Dashed lines are cor- 
responding predictions from reference 2. 



Figure 5. — Weibull plot showing three-point and disk specimen dy- 
namic fatigue strengths that were loaded at a rate of 1 MPa/s. Solid 
line corresponds to CARES/LIFE prediction for C = 0.82 (shear- sen- 
sitive criterion). Long dashed line corresponds to the disk pre- 
diction for C = 0.90_and the short dashed line corresponds to the 
disk prediction for C» 1 .0 (shear-insensitive criterion). 
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